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Abstract 

We study the competitive equilibrium of large random economies with linear activities using methods 
of statistical mechanics. We focus on economies with C commodities, N firms, each running a randomly 
drawn linear technology, and one consumer. We derive, in the limit N, G — » oo with n — N/C fixed, a 
complete description of the statistical properties of typical equilibria. We find two regimes, which in the 
limit of efficient technologies are separated by a phase transition, and argue that endogenous technological 
change drives the economy close to the critical point. 

1 Introduction 

The aggregation of microeconomic behavior into macroeconomic laws is a difficult task because of the pres- 
ence of heterogeneity both at the level of individual characteristics and of interactions. The paradigm of the 
representative agent, which essentially reduces the problem to that of a single macro individual, has shown 
all its inadequacy [1], calling for alternative approaches. Computational methods - both in the spirit of agent 
based modeling or implementing general equilibrium theory - represent a viable substitute, rapidly growing in 
popularity. However these techniques provide punctual results which are difficult to generalize. While they are 
very useful in deriving specific results for a specific economy, they do not lead to a broad understanding. 

At the other extreme, the methods of mathematical economics aim at general results - such as existence, 
uniqueness, efficiency - which hold for broad classes of situations. Pinning down the typical macroeconomic 
behavior beyond these general results is however very hard, especially when agents are heterogeneous (e.g. in 
their endowments, technologies, budgets, utility functions, . . . ) and are interconnected via a complex network 
of interactions. 

Understanding the complex macro-behavior of a system does not necessarily require a detailed description 
of it in all its complications. Indeed many laws which govern macro-behavior have a statistical origin. E.P. 
Wigner [2] first had the intuition that in such cases, the collective behavior of a large system with N degrees 
of freedom - heavy atoms in his case - is well approximated by that of a system with random interactions in 
the limit N — > oo. Indeed, if the relevant properties obey laws of large numbers, then they will be substantially 
independent of the specific realization of the interactions when N is large. 

The statistical properties of random systems have been a central research issue in statistical mechanics for 
the past two decades, and extremely powerful analytical tools to calculate them have been developed. These 
techniques have already found a wide range of applications outside physics: among others, in combinatorial 
optimization problems and computer science, in the theory of neural networks, in information theory, and in 
agent based models [3-7]. 

It has been realized several times by different authors [8,9] that tools developed in statistical mechanics can 
be useful in economic theory. While the idea of studying large random economies as a proxy for a complex 
economy with heterogeneous agents may not be entirely new (see e.g. [10]), modern tools of statistical mechanics 
of disordered systems have not yet been exploited. Here we apply these tools to the study of the typical properties 
of large random production economies. The model we shall consider, outlined in detail in the next section, is 
based on a C dimensional commodity space and has N firms with linear technologies, as in [11]. Feasible 
technologies are assumed to be drawn at random from some probability distribution. Firms chose technologies 



from the set of feasible ones and fix the scale of operations so as to maximize their respective profits. The total 
supply is matched to the demand of a single consumer with initial endowments drawn at random from a given 
distribution. Equilibrium prices, operation scales and consumption levels are determined by imposing that all 
markets clear. Equilibrium quantities are random, as they depend on the draw of technologies and o initial 
endowments. A complete statistical characterization of an ensemble of equilibria of large random economies is 
obtained in the limits N — > oo and C — > oo . The laws we derive are of a statistical nature, i.e. laws that a 
typical realization of an economy from the ensemble will satisfy almost surely, i.e. with a probability close to 
one when N is large. 

We show in particular that this approach: i) identifies the relevant macroscopic variables, the so-called order 
parameters, describing the behavior of the system in the limit N — > oo; ii) allows the calculation of the values of 
the order parameters from the solution of a "representative" firm problem, which embodies all the complexity 
of the full heterogeneous model; iii) enables one to derive distributions of consumption levels and of scales 
of activities at equilibrium. We will prove that for a broad class of choices the properties of the competitive 
equilibria change qualitatively at a critical value n c — 2 of the ratio n = N/C. This change becomes a sharp 
phase transition in the limiting case of efficient technologies. Loosely speaking, the economy expands rapidly 
when n increases for n < 2 whereas, when n > n c the economy is in a mature phase where the technology space 
is to a large extent saturated. Even though our picture is static, we shall claim, in the final section, that in a 
dynamic setting technological innovation (i.e. changes in N and/or C) driven by total output growth are likely 
to drive the economy close to the critical value N/C = 2. 

After discussing the model in the next Section we present, in Section 3, the main results. In order not 
to obscure the emergent picture, a detailed account of the approach and of the calculation is given in the 
appendix. More specific and quantitative results will be discussed in Section 4 and in Section 5 we argue that 
economies self-organize close to the critical point n ss 2. We close by summarizing our results and discussing 
some perspectives in the final section. 

We made an effort to keep the discussion and the mathematical complexity at the simplest level, even at 
the price of introducing restrictive or unrealistic assumptions. The present approach can however be easily 
generalized to more realistic (and more complicated) models. 



2 The model 

We consider an economy with C commodities, N firms endowed with random technologies, and one represen- 
tative consumer. Firms strive to maximize their respective profits, while the consumer aims at maximizing 
his utility. The two problems are interconnected by the market clearing condition. In detail, we consider the 
following setup. 

The company i (i = 1, . . . , N) is characterized by a technology (or activity) with constant returns to scale 
that, when run at scale Sj = 1 produces qf > or consumes qf < units of commodity c (c = 1, . . . , C). 
If the technology q i = {qf}c=i is operated at a scale s» > 0, then firm i produces or consumes Siqf units of 
commodity c. Technologies cannot be reversed, i.e. s, > 0. Following [11], we do not restrict our attention 
to Leontiev input-output models: each activity can have several outputs (joint production) and there are no 
primary production factors (i.e. qf > is possible for all c). As in [11], it will be important to impose that it 
is impossible to produce a positive amount of some commodity without consuming a positive amount of some 
other commodity. A sufficient condition to ensure this is that 

c 

Vz: $> c = - £ . (1) 

c=l 

Here e is the difference between the quantities of inputs and outputs, which measures the inefficiency of the 
transformation process of technology i. 

The profit of firm i is 7Tj = Sj(p • q i ) = Sj J2c=i QiP c > where p = {p c }^ =1 is the price vector, which we assume 
to be non-negative. Each firm fixes Si by solving the problem 

max TTi (2) 

Si >0 

at fixed prices. 

The representative consumer, whose utility function we denote by U(-) and whose initial endowment wc 
denote by x , chooses his consumption x = {x c }^ =1 by solving 

max U(x), B = {x c > : p ■ x < p ■ x } (3) 
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at fixed prices. B represents the set of consumption plans that satisfy the consumer's budget constraint. 

At equilibrium, the total supply of each commodity is required to match the demand from the representative 
consumer (market clearing), i.e. 

N 

Vc : x c = Xq + s ilt (4) 
i=l 

The simultaneous solution of the maximization problems (2) and (3) subject to (4) constitutes the competitive 
equilibrium we will study in this paper. 

Before specifying further our model, it is worth to make a couple of remarks. First, multiplying both sides 
of (4) by p c and summing over c one finds that, in equilibrium, 

N N 

p- (x - aj ) = ^Si(p- = ^7Ti = (5) 
i=i i=i 

The last equality comes from the fact that p ■ (x — x ) < because of the budget constraint and ni > because 
firms can always achieve iri = by not producing. So, on one side one recovers Walras' Law p ■ (x — xq) = 0, 
while on the other we find that iri = for all firms. 

Notice also that, combining (1) with the market clearing condition we find that in equilibrium 

c N 

5> c -a:g) = -e$> (6) 

C=l 1=1 

This equation means that total equilibrium consumption will be lower than the initial one. The model thus 
focuses on the ability of the productive sector to provide scarce goods (with small Xq) using as inputs abundant 
commodities (with large Xq) so as to increase welfare. 
We assume that technologies qf are given by 



1 C C 

C=L 

where are independent Gaussian random variables, with zero mean and variance A/C and the last two terms 
enforce the constraint (1). Appendix B shows that the assumption on the distribution of q\ can be relaxed 
considerably for our purposes 1 . 

In what follows, we shall use the notation (. . .) for expected values over the distributions of the 

variables u, v, . . . , z, but we shall omit the subscript when no confusion is possible. 

Commodities are a priori equivalent. The initial endowments xg are drawn at random from a distribution 
p(-), independently for each c. Furthermore we shall also suppose that 

c 

U(x) = "£<x c ) (8) 

c=l 

where u(-) is postulated to be increasing (u'(x) > 0) and convex (u"(x) < 0). These assumptions, which simplify 
our analysis considerably, appear to be extremely restrictive. They appear less unrealistic considering, as in [12], 
that x may measure desirable characteristics or properties of commodities rather than quantities thereof. In 
this light, the departure from Leontiev technologies with a single output becomes natural. 

It is also useful to introduce a measure of economic activity similar to the Gross Domestic Product (GDP). 
The total market value of all goods produced is the sum of (x c — xg)p c for all c with x c > Xq. Because of Walras' 
law (5), this is equal to \ J2 C \ xC ~ x o\p c - Normalizing prices to the average price level, we obtain 

GDP - c ^ 1 ^' X ^ pC ■ (9) 
2 E c =iP c 

The key parameters of the model are thus N, C, e, A, the distribution /?(•) of the initial endowments, and 
the utility function u(-). We shall focus on the non trivial limit N oo defined as 

(n) 

lim = lim (10) 

n=N/C 



1 In fact, any distribution which satisfies (1) and has a characteristic function log 

(eikqf^ = ^(k/VC) with iP(x) = ~Ax 2 /2 + 

0(x 3 ) would leave our results unchanged. 
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where n = N/C is held fixed as N — > oo. 

A simple geometric argument, for e = 0, suggests that n — 2 will play an important role. Let us write the 
initial endowments as Xq — x n + Sx^, separating a constant part (x ) from a fluctuating part (Sxq) such that 
J2 C Sx c = 0. With e = 0, Eq. (6) implies that the component of consumption along the constant vector remains 
constant. All the transformations take place in the space orthogonal to the constant vector: q i ■ x = q i ■ Sx . 
In other words, those technologies with q i ■ Sx < which reduce the initial spread of endowments Sx lead to 
a increase in wealth and hence will be run at a positive scale. Those with a positive component along 6xq will 
have Si = 0. Given that the probability to generate randomly a vector in the half-space {q : q ■ 8x < 0} is 1/2, 
when N is large we expect N/2 active firms. Still the number of possible active firms is bounded above by C, 
hence when n = N/C — 2 the space of technologies becomes complete and x c = Xo Vc. There is no possibility 
to increase welfare further. We shall see that n = 2 separates two distinct regimes of equilibria even with e > 0. 

It is easy to see that the problem of finding equilibrium prices, production scales and consumption levels of 
the economy in the above setting is reduced to the following: 

max U Ixo + ^ SiqA (11) 

Given the solution {s*} to this problem, the equilibrium consumption levels are given by the market clearing 
condition (4) and the (relative) prices are derived from marginal utilities as 2 

= u'(x* c ) (12) 

Eq. (11) is a typical problem in statistical mechanics. The general approach to this type of issues is discussed 
in Appendix A. Equilibrium quantities are random variables because of the randomness in the technologies q i 
and in the initial endowments xo- Still there are statistical properties of the equilibrium which hold almost 
surely in the limit (10). These will be the subject of our interest. In Section 3, we present the general solution, 
while in Section 4 we shall specialize to specific examples. The reader interested in technical details is referred 
to the appendices for a detailed account of the method and of the explicit calculation. 



P = 



dU 

dx c 



3 The solution and its generic properties 

As shown in appendices A and B, the solution of the equilibrium problem (11) in the limit N 
n = N/ C fixed is given by 



where 



W I I ( _ \\ 



with 



(13) 



h(fl,K,p,a,x,X) = ( max 

\ s>0 



(ta - ep)s - - X s 2 



-n X +-K P - — — - 
2 n 2nA 



2n 



XV + 



H — ( max 

n \ x>o 



(x) — — \ x — xq + k + V nAQtj 



(14) 



t,Xo 



and fl*, . . . ,x* are the saddle point values of the parameters, i.e. those which solves the system of equations 
|^ = 0, . . . , || = 0. The variables O, n,p, a, x, X are called order parameters in statistical physics 3 . They 
emerge from the analytic approach (see appendices A and B) as the key macroscopic variables which describe 
the collective behavior of the equilibria. 

In Eq. (14) t is a Gaussian r.v., with zero mean and unit variance and as usual (. . .) t , (. . .) t Xg stand for 
expectation values on t and on t and xo, respectively. The precise derivation of this result is described in 
Appendix B. 

The structure of h is reminiscent of the original problem. The first term on the r.h.s. can indeed be regarded 
as the profit maximization of a "representative" firm. The variable s is indeed one of the variables Sj which 

2 Utility maximization under the budget constraint, Eq. (3), yields dU/dx c = Xp c where A is the Lagrange multiplier imposing 
the budget constraint. We can take A = 1 exploiting the invariance p c — > ap c for any a > 0, thus fixing the level of absolute prices. 

3 In order to keep notation simple, we shall generally omit in what follows the asterisk on order parameters which denotes saddle 
point values. 
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appear in the original problem (11). The solution of the maximization problem in the first term of Eq. (14) is 
given by 

(ta-ep)/x if t>ep/ a 
if t < ep/cr 



s*(t) = 



(15) 



Since t is a random variable, s* is also a random variable and its probability density can be derived from that 
of t. The result is 



Q(s) = (1 - cj>)6(s) + 



X 



2ir<r 



0(s) exp 



(xs + epf 



2a 2 



-erfc 
2 



V2cr 



(16) 



where 0(s) = for s < and 0(s) = 1 for s > 0. The variable s is the scale of production of a (representative) 
firm, hence (16) yields the distribution of s, in the economy and <f> is the fraction of technologies that are active 
(i.e. such that s, > 0). 

Likewise, the last term on the r.h.s. of Eq. (14) is related to utility maximization with respect to a 
"representative" commodity. The variable x is indeed one of the variables x c which appear in the original 
problem. The solution of this problem is given by 



x*(t, xq) : x u '( x *) = x * ~ x a + K + VnAfit 



(17) 



which is always positive provided u'(x) — > oo for x — > 0. The probability density of x c in the economy can 
be derived from that of t and xo in the same way as above for the scale Sj of production. The conditional 
probability of x c given x% is computed in Appendix C. The result is 



P{x\x ) 



1 - xu"{x) 



exp 



(x — xo — x u '( x ) + K Y 
2nMl, 



(18) 



Hence the variable x — x u '( x ) has a Gaussian distribution with mean xo — k and variance nAfi. 

The two "representative" problems are coupled in a nontrivial way through the other terms in (14). 

The structure of the solution becomes more clear if we analyze the set of saddle point equations |^ = 
0, . . . , || = 0, with 6 — (fi, K,p, a, x, X.)- After some algebra (see appendix B), these can be cast in the 
following form: 



P 

X 

a 

X 

K 



= <«V)> M 



t,XQ 



(u'(x*))' 



t,xo 



(u>(x*)) 



t,xo 



nA 



= PX + ne(s*) t 



(19) 
(20) 

(21) 
(22) 

(23) 
(24) 



The first of these equations relates the parameter p to the average (relative) price because of (12), while the 
third one implies that a is a measure of price fluctuations 4 . 

Using these relations (see appendix D for details) one finds that at the saddle point 



h(Cl*,K*, P *,a*,x*,X*) = -Hx*)) 



t,xa 



(25) 



This is indeed what we expect looking at the original problem (11). Furthermore, taking the expected value of 
(17) and combining it with (19) and (24) yields X (u'(x*)) txo = X P = ( x *) t .x - ( x o) Xo + K = ( x *)t,x ~ ( x o) Xa + 
PX + ne (s*) r Then 

( x *) t> x = ( x o) xo -ne(s*) t (26) 

which is exactly Eq. (6). Finally, it is possible to show (see Appendix D) that Eq.s (19-24) also "contain" 
Walras' law in the form 

(u'(x*)(x* - x )) t xo = 0. (27) 

4 It is possible to derive an explicit analytic form of (22), (23) and (24) in terms of error functions. The present formulas are 
however more suited for the discussion which follows. 
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The dependence on A of the solution can be clarified by a rescaling argument: changing variables to p' = p, 
X 1 = x/A, a' = ct/VA, ^' = Afi, x' — X an d K ' — K one finds that the solution only depends on the parameter 
e' = e/A. Hence the behavior of the solution with respect to A is easily related to the dependence on e with 
A = 1. Notice that, a dependence on A remains after the change of variables in the distribution of Sj, Eq. (16). 
This means that production scales satisfy the scaling relation 

Si (A) = Si (l)/\/A (28) 

The behavior of the solution when the spread of the initial endowments (Sxfy = ((xo — (xo)) 2 ) is very small can 
be computed with asymptotic expansion methods. The key observation in the expansion (see appendix E) is 
that x* also has very s mall flu ctuations. This, in turn, implies that prices also have very small fluctuations, 
indeed a = \u" ((xo))\-\/ (Sx'q) . The scales of production also vanish when (Sxq) — > 0, but with a singular 
exponential behavior: 

not(6x 2 Q ) V2 e- A /( Sx °), (29) 

for some constant A. Hence we find that no economic activity takes place (<p — > 0, — > 0, (s*) t — > 0) in the 
limit of uniform endowments (see Appendix E for technical details). This is what one should expect from the 
beginning: when the consumer is endowed with the same amount of equally valued commodities, there is no 
transformation (with e > 0) which can increase welfare. 

A further interesting limit, for which we can derive generic results, is that of vanishing e. Setting e = one 
finds in a straightforward way that £1 = er 2 /(2x 2 ), x — n ^/(^X) an d k = px- Eq. (7) yields <j> — 1/2 which 
means that half of the firms are active, in agreement with the geometric argument of the previous section for 
n < 2. When n — > 2~ the equations develop a singularity: Indeed x cx (2 — n), a cx \/2 — n vanish whereas 
the average scale of production diverges (s*) cx 1/^/2 — n. A detailed account is given in Appendix F. The 
case n > 2 is more subtle as it requires a careful asymptotic study of the limit e — > where again realizing 
that x* has small fluctuations of order e is crucial. The bottom line is that (see appendix F for details) price 
fluctuations vanish linearly with e, i.e. a cx e but also x cx e so the factors ep/x and <r/x in Eq. (15) are finite. 
Hence scales of production remain finite, as e — > and they diverge when n — > 2 + as (s*) cx — 2. The 

fraction of active firms turns out to be <p — 1/n, which means that there are exactly C firms operating. 

Eq.s (25), (26) and (27) show that the saddle point equations, which represent the simplest mathematical 
description of the random economy in its full complexity, manage to capture in a compact, though somewhat 
intricate, way the basic properties of the economy. This is a useful consistency check. The best way to unravel 
the resulting behavior beyond these generic laws is however to specialize to particular cases. 



4 The solution: typical cases 

In this section we display the behavior of the solution outlined in the previous section for some particular choices 
of the functions u(x) and p(xo). In spite of their apparent complexity, Eq.s (19-24) can be solved numerically 
to any desired degree of accuracy. Using the scaling argument above, we can safely restrict ourselves to study 
the dependence on e setting A = 1, without any loss of generality. 
We shall henceforth set 

u(x) = logx. (30) 

We start our discussion from the case 

p(x Q )=e- X0 , x Q >0. (31) 

Fig. 1 compares the numerical solution with computer experiments. We generate many realizations of the 
random economy and compute numerically the equilibria for each of them. The analytical results we obtain 
in the limit C — > oo turn out to give a quite accurate description of the behavior of relatively small systems 5 
(i.e. C = 16) even for a single realization. Fig. 1 shows that there are essentially two different regimes. For 
n < n c = 2 roughly half of the firms are active, whereas for n ^> n c the number of active firms saturates to C. 
The GDP also shows a similar behavior. It increases with n and saturates for n > 2. 

The transition between the two regimes becomes sharper when e decreases and it gives rise to a singularity 
in the limit e — > 0, as we have seen in the previous section. This is clearly visible in Fig. 2, where we plot the 
behavior of various quantities as a function of n for different values of e. 

For n < n c the average scale of production (s*) increases with n. This means that, in this region, existing 
firms benefit from the entry of a new technology (i.e. if N — * N + 1, see later). This positive complementarity 
arises because the new firm increases the availability of inputs to other firms. 

5 We resorted to a simple iterative scheme to converge to the equilibria. This fails to converge properly for C or N too large or 
for e < 1. 
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4 6 

n=N/C 



Figure 1: Comparison between the analytic solution for N — > oo (full line) and equilibria of random economies 
computed numerically for C = 16 and 32. The parameters are e = 0.05 and A = 1 whereas initial endowments 
are drawn from an exponential distribution. Dots refer to a single realization with C = 32 whereas the dotted 
(dashed) line is the average over 100 realizations for C = 16 (32). Left: n<p, which is the number of active firms 
(sj > 0) divided by C, versus n. Right: GDP versus n. Insets in these figures show the behavior of <j) and of 
GDP for e = 0.01, 0.05, 0.1 and 0.5 from top to bottom. 




Figure 2: Behavior of equilibrium quantities as a function of n for e = 0.5, 0.1, 0.05 and 0.01. In all cases 
A = 1 and p[x) — e~ x . Top left: (s*). Bottom left: relative price fluctuations. Here we identify prices with 
marginal utility p — u'(x*) and Sp = u'(x*) — (u'(x*)). Top right: average consumption (x*). Bottom right: 
relative fluctuations of consumption Sx — x* — (x*). 
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Figure 3: Probability density functions of operation scales s (left) and of consumptions x (right) at equilibrium 
for e = 0.01 and n = 0.5, 2 and 5. 



For n > n c , instead (s*) decreases with n, the introduction of a new technology typically causes a reduction in 
the scale of activity of the already existing firms. When e — > the curves develop a singularity (s*) ~ 1/ yfn — 2| 
at n c = 2, as discussed in the previous section. 

As n increases relative price fluctuations decrease. But the decrease becomes very sharp close to n c for 
e <C 1. In this case, at n c price fluctuations suddenly drop to a level close to zero. This is related to the 
behavior of the variable x* shown in the right panel of Fig. 2. In the region below n c the average consumption 
level decreases. In this region firms take advantage of the spread 



(Sx) __ V^ 2 )-(z) 2 ( . 

~W ¥) (32) 

between scarce and abundant goods to make a living. But as n approaches n c , the spread in x quickly drops to 
a very low value, making life more difficult. Increasing n beyond n c = 2, the economy becomes very selective 
toward increasingly efficient technologies which can perform the desired transformation between commodities 
with a smaller decrease in the average level (x) of consumption. This is clearly shown in Fig. 3 where we plot 
the probability densities of x for three different values of n = 0.5, 2 and 5. The right plot shows that while for 
n = 0.5 the distribution P(x\-), Eq. (18), retains the character of the distribution of initial endowments p(xo), 
it becomes more and more peaked around (xq) as n increases. At the same time the distribution of s, Eq. (16), 
becomes broader and broader. 

The distribution of Sj, Eq. (16), which may be considered as a proxy for firm sizes, gets broader and broader 
as n increases. Interestingly mature economies, such as Japan [13] or the US [14], are characterized by a very 
broad distribution of firm sizes, which we can put in relation with Q(s). The shape of the distribution found 
empirically is close to a power law, which is different from (16). However, it is not difficult to derive a power 
law distribution of s relaxing the unrealistic assumption that all firms have the same value of A and e [15]. 

The generic picture of the overall economy depicted thus far remains unchanged for different distributions 
p(xo) of initial endowments or for different utility functions u(x). For example, Fig. 4 shows the results obtained 
with 

p(x ) = (l-f)5(x )+fS(x -l). (33) 

This captures the situation where only a fraction / of the commodities is present in initial endowments (primary 
goods) whereas the remaining commodities have to be provided by the productive sector. The behavior of cj>, 
(x*) and relative prices is very similar to that found for the previous model. Fig. 4 shows that the average scale 
of production and the relative fluctuations of x c show a qualitatively different behavior. Again the two regimes 
with clearly distinct properties can be identified for n < n c and n > n c . 
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Figure 4: Left: Scale of production (top) and consumption fluctuations (bottom) as a function of n for a bimodal 
distribution of initial endowments (Eq. 33) with / = 0.2 and xo — 5 (e = 0.01). Right: distribution of x for 
n = 1, 2 and 5 from top to bottom. 



5 Discussion 

The behavior of the solution with n allows us to identify two classes of economies: mature economies (n > 2) 
with a full-blown repertoire of technologies which closely saturate consumer's demand and immature economies 
(n < 2) characterized by few technologies scattered in a large space of productive opportunities. 

Strictly speaking, our considerations must be limited to comparative statics in view of the static nature of the 
equilibrium we study. However it is suggestive consider dynamic transitions between equilibria. In particular, 
a transitions N — > N + 1, corresponding to the draw of a new technology, can be considered as the result of 
the discovery of a new method of combining inputs to produce desirable outputs; a new design as Romer [21] 
calls it. Note however that in Romer's model innovation entails the discovery of a new intermediate commodity 
and there is no real "combination" of inputs and no heterogeneity across technologies. By contrast, innovation 
in our model describes the expansion of the frontier of feasible industrial transformation processes by discovery 
of a new activities which is structurally different from existing ones 6 ; both the discrete nature of designs and 
the uncertainty of the discovery process are retained. Whether an innovation leading to the draw of a new 
technology is adopted or not will depends on the specific technologies that are already present 7 Generalizing, 
one may also consider transitions C — ► C + 1 introducing a new commodity or characteristic [11], or which 
makes it possible to consume it. 

In such a dynamic view 8 , it is essential to consider the incentives for innovations in order to understand 
which transitions will most likely be generated endogenously. If we assume that transition rates depend on 
investments in research for new technologies, and that investment are related to expected changes in GDP 
generated by technological changes, then Fig. 5 suggests that the economy will drift towards n « 2. Indeed 
transitions N — > TV + 1 cause an increase in GDP which is sizeable for n < 2 and almost negligible for n > 2 
(specially for small e). On the contrary, transitions C — > C + 1 which decrease n increase substantially GDP 
only for n > 2. 

The same conclusion can be reached assuming that investment in research arise from the productive sector 

6 This is only one of the possible modes of technological innovations. Innovations may also increase the efficiency of an existing 
technology, e.g. decreasing the input requirements, which may be captured by changes in ej and A;. Our focus here is on structural 
technological change. 

7 It has been argued that technological innovation is path dependent [16] . This means that the draw of the N + 1 st technology 
depends on the N technologies which are already present. Such issues can clearly not be addressed within our quasi-static approach. 

8 For example consider the following cartoon of an evolving economy: In each period t Nature endows the representative consumer 
with a bundle xo(t) of commodities, drawn at random from the distribution po- Given the existing technologies, {q i : i = 
1, . . . , N(t)} this is transformed into the optimal bundle x(t) by the productive sector then x(t) is consumed at the end of period 
t. Finally an innovation event, i.e. a transition N(t) — > N(t + 1) = N(t) + 1 or C'(t) — » C'(t + 1) = C(t) + 1 may take place. 
It is implicit in this description that technological change occurs on a much longer time scale than that needed for the economy 
to reach equilibrium. This view may also describe the way in which existing technologies diffuse in a developing country. Then 
technologies are not discovered anew, but just become operative or feasible as e.g. human capital or infrastructure accumulates, or 
as institutional constraints are removed. 
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Figure 5: Variation of GDP for N — > N + 1 (full lines) and for C — > C + 1 (dashed lines). Here u(x) — log a; 
and p(xo) — e~ xa . 



itself. Indeed the average scale (s*) of activity of firms increases with n for n < 2, which means that a transition 
N — > AT + 1 causes an increase in the average scale of activities already active. This means that, at the 
equilibrium prices of the economy with iV < 2C technologies, the profits of already existing firms increase on 
average when the new technology is introduced. Likewise, the decrease of (s*) with n for n > 2 suggests that 
transitions C — > C + 1 increase firms profits, on average, for n > 2. This again yields a drift towards n ~ 2 due 
to endogenous technological change. 

These arguments suggest that economies may spontaneously evolve towards a critical state, i.e. that they 
may be a further realization of self-organized criticality [17]. 

6 Conclusions 

Summarizing, we have addressed the problem of calculating the general equilibria of large linear production 
economies with random technologies and a single consumer with tools of statistical physics. In a nutshell, 
our results can be stated as follows. When the ratio n of the number of available technologies to the number 
of commodities is below a threshold n c = 2, the average operation scale grows as n increases and roughly 
one-half of the firms are active. For n < 2, new technologies are easily accepted and the economy on the 
whole expands with n. When n > n c , instead, the production sector is saturated, i.e. the number of active 
technologies converges to the number of commodities, and new technologies are accepted only at the expense 
of reducing the operation scales of the other technologies. The transition becomes more and more sharp as the 
parameter e, measuring the inefficiency of each technology, approaches zero. From the consumer's viewpoint, 
welfare increases with n in both regimes. The main component of welfare increase with n is different in the 
two regimes: for n < n c , welfare level grows with n because the spread in consumption levels (Sx) decreases 
with the introduction of technologies which transform abundant commodities into scarce ones. For n > n c , 
instead, growth arises from the introduction of more efficient technologies, granting an increase in the level of 
consumption (x*). Accordingly, the relative spread of prices (Sp) / (p) decreases with n. 

Considering the incentives for technological innovations, we uncover a mechanism by which the economy 
self-organizes to the critical state w 2. Our model clearly is unrealistic in many respects. Still it may capture 
some novel aspects of structural technological change. The extension of these approach to a fully dynamic 
setting, including capital accumulation, may shed new light on theories of endogenous economic growth. 

Above all, we propose the use of statistical mechanics of disordered systems to study the typical properties of 
the general equilibria of large random economies. We have shown how these methods are able to deal effectively 
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with heterogeneity, providing a complete statistical description of the equilibria, which is consistent with generic 
results. The relevant quantities - called order parameters - are naturally identified by the method. Given the non 
standard type of the calculations involved, we also present computer experiments which convincingly support 
our results. 

The approach generalizes in a straightforward way to more complex situations and we hope that this work 
will stimulate cross-fertilization between the fields of economic theory and statistical mechanics. 
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Appendices 

A The method 



The standard technique to maximize a function of N variables with TV — > oo in statistical mechanics relies on 
the well-known steepest descent, or saddle point, method. Let Hn(-) be an extensive 9 function of s = {si}fL l7 
and imagine that we want to compute the maximum value of Hn(s)/N in the limit N — > oo. Then 



lim — max Hm(s) = lim lim 

W^oo N s /3-.oo N 



1 



~w XogZN 



where 



Z N {/3) = J ds 



(34) 



(35) 



is called the partition function associated to H jy. Here J ds stands for an iV-dimcnsional integral on the whole 
domain of definition of s. The idea of (35) is that the integral for (3 3> 1 is dominated by regions where Hm 
is maximal. This recipe turns the problem of maximizing h into that of calculating Zm and evaluating the 
asymptotic behavior of its logarithm. 

This task becomes much more difficult when Hm depends on a set of random variables q with probability 
density p{q). We denote this dependence by H^{-\q). The generic situation is that q enters the definition of 
the interactions among the N components of s and Hm is a sum over all interaction terms. In such situations, 
we expect that a sufficiently regular Hm will obey the law of large numbers, so that e.g 



lim — max Hiy(s\q) = lim — 



max Hjs;{s\q) 



(36) 



In other words, max Hm/N is expected to be a self-averaging quantity, namely to have vanishing sample-to- 
sample fluctuations in the limit N — > oo. If one wanted to generalize (34) to the evaluation of (36), one would 
have to compute the g-average of the logarithm of the partition function Zjy(l3\q)- Unfortunately, the logarithm 
prevents every useful factorization of such an average and makes this way impracticable. 

The replica method is the standard statistical mechanical technique to circumvent this difficulty. Using the 
formula 

log Z N ((3\q) = hm 

i — >o r 

we can reduce our problem to that of computing ([i?Ar(/3|q)] r ) . 
it amounts to computing 



(37) 

f 

This is feasible for integer values of r because 



[z N (p\q)y 



e PH N (s\ q ) ds 



1 H N (s a \q) 



n ds 



(38) 



9 I.e. such that there are two constants k+ and fc_ satisfying k-N < Hm(s) < k+N for all N and s. 
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which is the partition function of r "replicas" of the original system with the same disorder realization q (hence 
the name of the method). The last step consists in performing an analytic continuation for real values of r and 
taking the limit r — > 0: 

lim 1 (max H N {s\q)) = lim lim lim l -L log {[Z N [fi\q)] r ) (39) 

The existence and uniqueness of the limit r — > 0, which looks somewhat bizarre, have been much debated in the 
physics literature (see [18] for a discussion). Even if this method remains a formally non-rigorous procedure, 
several rigorous mathematical results confirm its validity in problems that are more complex that the one we 
deal with here [19,20]. We hope this (together with the agreement with computer experiments) gives the reader 
a sufficient level of confidence to accept the r — > passage. 

The technical part of the calculation lies in the introduction of a finite number of auxiliary integration 
variables = {9\, . . . , 9k} allowing the averaged replicated partition function to be re-cast in the form 



(lZ N (p\q)Y) q - J e? Nr l h W + °MWd0 



(40) 



where h(-) is some function and o(r, (3, ■) — * in the limits (3 — > oo, r — > 0. The 9 variables are called order 
parameters. Their nature and number are dictated by the mathematical structure of the problem (see Appendix 
B for the details of our case). Finally, assuming that the limits r — > and N — > oo commute, the latter can be 
taken first in (39) thus allowing to evaluate (40) by the saddle-point method as 

([Z N (J3\q)] r ) q ~ e P"r[h{6*)+oW)] (41) 
where 0* is the saddle point value of 9 which dominates the integral in (40). Hence, putting things together, 
lim 1/max H N (s\q)) = lim lim lim — L log e pNr ^ ' = h t 9 *^ ( 42 ) 

The core of the procedure lies in (40) where, by a lengthy calculation one identifies the relevant order 
parameters and the function h. This crucial but technical step is presented below (appendix B) for our 
problem. 



B The explicit calculation of the representative agent problem 

The partition function in our case reads 

Z N {(3\q) = [°° e 0u ( Xo+ ^ Siq ^ds 
Jo 



(43) 



with U(x) the utility function of the representative consumer. As stated above, in order to analyze the statistical 
properties of the equilibria we have to evaluate {[Zjs[(f3\q)] r ) q and resort to (39), with Hn given in our case 
by U and with all the necessary constraints. Before proceeding, we shall introduce some useful definitions and 
identities. The first one is the (5-function 5(x) which is defined through the relation 

f(y) = [ 6(x - y)f(x)dx (44) 

JR 

for any function /(•) and y <G R. We will also use the exponential representation of the (5-function 



Another mathematical tool we will use is the Gaussian or Hubbard-Stratonovich transformation, viz. 

dx 



exp 



"6 2 " 




\ x 2 , 1 


2 


= / exp 

JR 


- Y + bx 



2ir 



(46) 



which allows to linearize arguments of exponentials at the cost of introducing averages over Gaussian r.v.'s. 
Now, in order to perform our calculation it is convenient to replace the consumption variables x by writing 
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explicitly the market clearing condition (4) in the partition function (43). To do so we use the defining property 
(44) of ^-distributions and write 



/>oo />oo C N 

Z N ((3\q) = dx ds e^W J] 5[x c - x% - £ s t qt) 

^° ^° c=l i=l 



(47) 



As already explained in Sec. A we will have to take the following steps: (a) average the partition function of r 
replicas over technologies, as in (40); (b) identify the correct order parameters of the problem to write the latter 
average as in (40); (c) take the limits N — > oo and r — > and get something of the form of (42); and finally (d) 
find the values of the order parameters at the competitive equilibrium (i.e. when j3 — > oo). 
The partition function of r replicas reads 



[z N (p\q)\ r = / ]ldx a / Uds a 

J a =l J a =l 



a— 1 c— 1 



N 



(48) 



Notice that the dependence on the technologies appears in the market clearing condition only, so that the 
average (. . .) involves only the last part of [ZN((3\q)] r . This average one must take into account the constraint 
(1), i.e. 



<-)« 



nr=i<K£c=i<tf+e) v 



(49) 



where (• • •)' stands for the average over unconstrained i.i.d. Gaussian vectors q with zero mean and variance 
(q 2 ) q = J2 C ((l C ) 2 ) q = Using Eq. (45) for the constraints, the denominator becomes 



N 



\i=l c=l / , 

while for the numerator we get 

INC r C N * 

U S (zZli+ e )IiIi S { xC a- xC 0-zZ 



N 1 



exp 



2 -, 



e 
2A 



(50) 



\»=1 c=l 



i=l 



r N d? r r d9 r N r ° A N C / r \ 

I n £ I jj d^ ^ j k E % + . E E ^ ^ _ ^ _ A_ E E ^ _ E ^ 

-i— 1 a— 1 i= 1 a— 1 c— 1 i— 1 c— 1 a— 1 



(51) 



Notice that the expected values involved in these calculations are all of the form ip(y) = (e %vqi ) . This is the 

characteristic function of qf and for the assumed Gaussian distribution it takes the form ip(y) = e~ Ay2 /( 2C ). 
This result can however be extended to any distribution of qf with ip(y) = ip(y/^/C) with a leading behavior 
ip{x) = — Ax 2 /2 + 0(x 3 ). Indeed all higher order terms in the power expansion of ip give vanishingly small 
contributions with respect to the first, in the limit C — > oo. 
Gathering all the terms, we have 



i=l Z7T J o=l Z7T J ° a=l J =1 o=l i=l 



r C 



N C 



(7 = 1 C=l 



i=l c=l 



" 1 

TT 1 



exp 



e 

2A 



(52) 
(53) 



In order to write the above in a form as simple as (40), the set of order parameters to be introduced must 
allow for a decoupling of the integrals over the variables %, Si^ a and x c a in such a way that the integrals on the 
different variables can be factorized. Here it is enough to introduce the following order parameters 



1 N 

'ab ^ ^ s i,a s i,b 



and 



1 N ^ 



i=l 



(54) 
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through identities such as 



1 = / du ab NS (N^ ab = J ^^e^N^-E^ 



iSi,b] 



(55) 



Then last term in the exponent of the numerator of Eq. (53) becomes 



N r 2 N r r 

E(*-E =E^- 27V E k ^a + N J2 



(56) 



i=l 



a=l 



a=l 



a, 6=1 



This allows us to separate the problem into three parts. Indeed we can re-cast the replicated partition function 
in the form of a set of integrals over the order parameters: 



([ZN((J\q)] r )< 



n 



dlOabdWab 



IT ot°M° exp [ Nh iWab}, {^ab}, {k a } , {k a })] 



11 Airi/N J 11 2ml N 

a, 6=1 ' a=l ' 

where h — g\ + .92 + .93 is the sum of three terms which can be computed independently. In particular 

.91 = fll({Wa&},{£ o6 },{fc a },{fc a }) = -9 E ^abUab-^2k a k a 



(57) 



(58) 



,6=1 



1=1 



.92 = .92({w Q {,},{fc a }) = log J 7^ j Yl dSaCxp [^ E ^a6S a Sfc + 2^fc a s a + iez- ^ 

-log 



A. 



0=1 



0=1 



0,6=1 



exp 



e 

2A 



(59) 



1 /" dx r°° ' 

.93 = .93({^a6},{fca}) = ^E l0 S / II / II ^ ^ E + 1 E ^ (*» ~ X o) 

c=l a=l ^ ^° a=l a=l a=l 

r r 

E X a X b UJ ab +nA^2 ^afcaj 



nA 
"~2~ 



(60) 



a, 6=1 



with n = N/C. The order parameters fc a have appeared after using an identity similar to (55) for k a . Now (57) 
is precisely of the form (40). 

In the limit N — > oo the integrals appearing in (57) are dominated by the contributions coming from the 
saddle-point of h and the solution of our specific problem can be written as 



(") 1 / \ 1 

Jim / m ax U(x)) = lim lim l— h({u)* b }, {u)* ab }, {k*}, {k*}) 

AT— >oo iV \ S lq P^ool — >0 fjr 



(61) 



where the * means that parameters take their saddle point value, i.e. those who solve the system of equations 

dh 



duj a l 



dh dh n dh 

= 0, 7^=0, ttt-^0, -^=0 

duj ab dk a Qka 



(62) 



for all a, b = 1, . . . , r. Ideally one should first solve these equations for generic r and then take the limit r — > 0. 

A word about the meaning of the order parameters introduced thus far is in order before taking the limit 
r — > 0. Indeed, u> ab is an r x r matrix for integer r, but it is not clear how can we handle it in the limit r — > 0. 
When we replicated the partition function passing from (43) to (48) we essentially passed from a problem in 
which U(x a ) is to be maximized to an equivalent problem in which ^2 a U(x a ) is to be maximized. The latter 
sum is evidently left unchanged by a permutation of the replica indexes 1, . . . , r. Hence it must be expected 
that, as long as there is a unique maximum (as in this case), replica permutation symmetry is preserved also 
by the solution of Eq.s (62). Then we expect a solution of the form 



w a6 = Mtiab + W(l - Sab) 
U*ab = ^S ab +Uj(l - S ab ) 

k n 1=1 k k r — k 



(63) 
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This is the so-called replica-symmetric Ansatz, which simply expresses the conservation of the permutation 
symmetry. When multiple maxima with different statistical properties exist, this Ansatz fails because replicas 
can converge to maxima with different properties, and hence replicas are no more equivalent. This situation is 
ruled out in our case by the nature of the function we want to maximize. 

With Eq.s (63), it is easy to find an analytic expression of the function gi, g 2 and gz in terms of r and to 
perform the limit r — > 0. Substituting (63) into the definitions of g\, g 2 and #3, after some straightforward 
algebraic manipulations one finds 



lim -g 1 

7 — >o r 

lim -g 2 = 

r^O r 



i ^00 — LOtO^j 



kk 



lim -g 3 

r^O r 



log 
log 



ds e? v ^ 



dx eW^.-o) 



t,Xo 



where 



0W(x;t,x o ) = 0u(x) 



0V(s;t) 



(x — xq + \fnKuit — \nAk^ ^ 



2nA(n - to) 



- log[27rnA(fi 



-s* + 



1/2 T e 



(64) 
(65) 
(66) 

(67) 
(68) 



We must finally evaluate the limit — > 00. In this limit, a somewhat special role is played by the quantity 
X = - to). Notice that 

1 N 

tt-u=—^2(s t ,a-s l , b ) 2 (69) 

i=l 

is the distance between two replicas. The two vectors s a and Sb both converge to the unique solution of the 
maximization problem as — > 00. Hence, we also expect the distance ft — to to vanish in this limit. But looking 
e.g. at W/0 one realizes that in order to avoid annoying divergences or trivial limits this quantity must vanish 
in such a way that the product /?(ft — to) stays finite. In other terms, one wants that ft — to <~ 1/0 for large 0. 
If this is the case, the maximization problem has a well-defined solution. Hence we assume that lim^j — > oc x is 
finite. Similar arguments lead to the introduction of the following re-defined order parameters, which remain 
finite as — > 00: 



X = nA0(Q-to), 



ft 



X 







-inAfc 



(70) 



[ A0' ^ = r2Q 



(71) 



Inserting these in the previous expressions we find that the r.h.s. of (61) (which we for simplicity denote again 
by h) can be written as 



1 / 7x \ 1 1 / r 



ds 



n0 



log [°° dx e^te*,** A (72) 

JO I t,x 



where now the functions V and W road 



W(x; t, x ) = u(x) 



V(s;t) = -^ 



(x — x + k + ^JnAVLty 
2x" 



(tVl ~ A/? 2 + Ke)s 



(73) 
(74) 



We neglected the last term in W because it is vanishingly small in the limit — > 00 when \ is finite. 

When — > 00, again by steepest descent reasoning, only the maxima of V and W contribute to the integrals 
over s and x. Therefore we can write the final expression for h as 



h(Q,K,K,i,x,x) = (max 

\ s>0 




(x — x a + k + y/nAQty 



u(x) 



2x 



(75) 



t,Xo 
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The difference between this expression and the one appearing in (14) is again a trivial redefinition of the order 
parameters. If we let now x*(t, Xo) and s*(t) be the values maximizing the functions W and V, respectively, 
and therefore given by (15) and (17), we can then expand (75) to obtain 

h(n, k,k,j, x ,x) = -f <(**) 2 ) t + Vi-^ 2 (ts*) t + Ke (s*) t + ^) - h K + 

The last step is to derive the saddle-point equations from which the values that the order parameters take on 
at equilibrium can be calculated. Computing the derivatives of h with respect to the order parameters we get 



^T\ ^(( x * ~ Xo + K + V^Mlt)t) (77) 



dh 

dQ 2 A 2xVnO\ v *" ~ u 1 1 * '" '"/ t ,xo 

dh 



——k — (x* - x + k + VnAilt) 



dn n nx \ It 



XO 



(78) 



S - ; ^ (ts*) t + e(s*) t -- K (79) 

I - i^m'^-^K < 8 »> 

I = -^ + ^((**- i -»+< t +^ Tt > 2 L < 81 » 

| = -i«^>, + in (82) 



Using the relation (17) and setting p = —k, g — y/j — An 2 we finally arrive at Eq.s (19-24). 

C The p.d.f.'s of s and x 

We illustrate here the procedure for calculating the conditional probability density of x (the equilibrium con- 
sumption) given x (the initial endowment). The derivation of the distribution of s follows exactly the same 
lines. One can start from the identity 

/OC J i 

e- t2 / 2 6[x-x*(t,x )} (83) 
-oo V27T 

Then one can make use of the property 5(x — x*) = \f'(x*)\5[f(x)], where f(x) is a function with an unique 
root in x* . From (17), we take 



so that 



P{x\x )= f 

J — ( 



dt -t 2 /2 l ~X u "i x ) 



x — xq — x u '( x ) + k 

e " ' ' o 4 



V2^ VnAQ, 



t + 



VnMl 



(85) 



From this, taking the integral over t, one immediately finds (18). 



D Calculation of h at the saddle point and derivation of Walras' law 

Substituting s with s*(t) (15) and x with x*(t, xq) (17) we can re-write h as 

+ - « x *))t xn - 7T~ (i x * ~ x ° + k + ^ntf) (86) 
n t,x ° 2nx \ lt,x 
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Now it's a simple algebraic problem. For the first term on the r.h.s. we use (23); for the second and the fifth 
we use (24); the third and fourth cancel because of (22); finally for the last term we use (17) and then (21) to 
find, finally, 



2n X V ' lt.x Q 2nV v " / t,x 2n\A F 



(87) 



(25) follows immediately. 

In order to derive Walras' law, we note that when computing (s*t) v one can make the substitution t = 
(xs* + ep)/a (which is only valid when s* > 0). Then (23) becomes 



X = ^(nxQ + en (s*) t p) = (n\Q. - p 2 X + hp) 



(88) 



where we have used (24) in the last equality. Likewise, we can substitute t in the average of (20) by solving 
(17) for t. This yields 

n X n = x({u\x*)f) - k p-(u'(x*)(x-x )) tiXo (89) 

V / t,xo 

which can be substituted back in (88). This yields the desired result (27). 

E Almost uniform initial endowments 

In this section we study the limiting behavior of the economy when the spread of initial endowments is van- 
ishingly small. In particular we show that when the initial distribution of endowments becomes uniform the 
volume of productive activity vanishes. We take A = 1 for simplicity. We take xo = xq + Sxq , with xq a fixed 
value and 5xo a small random variable, and discuss the solution to the leading order in (<5x 2 ). Then, taking 
x* = xq + Sx* we can write, to leading order in Sxq and Sx* 



Xu'(x*) = x u '{x ) + xu"(xq)Sx* = x* — x + k + VnTit = Sx* — Sx a + k - 
From here we can identity the zero and first order terms in Sx* , viz. 

Sxq — VnTit 



'nClt 



K = X u '( x o), Sx* 



1 - xu"(x ) 



Then from Eq. (20) we get 
and from Eq. (21) 



X 



-U"(XQ) 

1 - Xu"{xq) 



u"(xq) 



1 2 



1 - xu"(x ) 



{(Sx f) xa + nO\ = x 2 [((^o) 2 )^ + nQ 



Coming to the equations for f2 and x we observe that s*(t) = s (t — r)9(t — t) where 



(90) 

(91) 
(92) 
(93) 



so = 



(Sx 2 ) 



Then we can write f2 = SqI(t) with I(t) = ((t — r) 2 9(t — t)) 4 , which can be solved for 



n = 



1 - nl(r) 



(94) 
(95) 



(96) 



Careful asymptotic analysis implies that in the limit of vanishing fluctuations (Sxq) — ► of the initial endowment 
t diverges as 

r = (97) 

u"{xq){5xI) 1 J 2 
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Using asymptotic expansion for 7(r) one finds that the leading order behavior of f2 is 

s i KM / (feo )2\ 3 / 2 e - 2[ «»w(i), . (98) 

Analogously for x we have 

/T , e 2 [V(- )] 2 

With these we finally get \ — \ u "(^o)\ and ° — \ u " (%o)\\J {$ x n) XQ - Therefore when the fluctuations of the initial 
endowment vanish there is no market activity. 



F Limit e -> 

Setting e = the averages over s*(t) become trivial and Eq.s (22,23) and (24) are easily evaluated. Progress 
with the other equations is possible, for generic u(x) and p(x a ), close to n = 2. Then we expect that the 

consumption vector x is nearly constant. Then, as in the previous case, we assume x* = (xq) + Sx* with Sx* 

small. Then p — u'((x )) as before. Using the expressions for O, \ and K an d expanding Eq. (17) to linear 
order as above, we find the expression 

5x* = (100) 

KI + & 

where u" = u"((xq)). This allows us to compute (u'(x*)t) = u" {tSx*) and hence to evaluate Eq. (20): 

X-AK|(l-|), n^2- (101) 

Likewise we can evaluate Eq. (21) and find a = \u"\W A(l — n/2). Using these in Eq.s (15,22) we find the 
divergence of = 1/(2 — n). 

This solution breaks down for n > 2. We need to take the limit e — > carefully into account. Again we 
anticipate that the spread dx* will be small. More precisely we assume that a and x vanish linearly with e and 
set so — cr/x and a — de and look for a solution with finite so and d. The existence of such a solution justifies 
our assumption. Then Q = Sg/(p/d) where 7(r) = ((t — r) 2 Q(t — t)) t has already been introduced above. The 
equation for x yields x — nAs n J(p/d)/(de) where J(r) = ierfc(r/v / 2). Expanding Eq. (17), as in the previous 
section and using the expressions for fl and x just derived, we find 



r * A x ° ~ ( x o))/ s o - y/nAI(p/d)t 

6X =6d nA\U»\J( P /d) (102) 

which justifies our a priori assumption of small fluctuations. Inserting this and the expressions of ft in Eq. (20) 
for x, after some manipulations, one finds 

nJ(p/d)leric(-^)=l (103) 

which gives dasa function of n and (x ). Notice that <f> — J(p/d) — 1/n. Furthermore J(r) < 1/2 for r > 
which means that this solution describes the region n > 2. This equation also simplifies the expression of 
X = As / (de) . The equation for a finally gives 



•>-{w£km (104> 

Note that Eq. (103) implies that d — > oo as n — > 2+. The leading behavior is d = j?^- In the same limit 

I{p/d) — > 1/2 which means that (s*) oc s ~ l/(n- 2) also diverges as n — > 2+ matching the divergence for 
n -» 2". 
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